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Abstract. We present the luminosity function (LF) of star clusters in M51 based on HST/ACS observations taken 
as part of the Hubble Heritage project. The clusters are selected based on their size and with the resulting 5 990 
clusters we present one of the largest cluster samples of a single galaxy. We find that the LF can be approximated 
with a double power-law distribution with a break around My = —8.9. On the bright side the index of the 
power-law distribution is steeper (a = 2.75) than on the faint-side (a = 1.93), similar to what was found earlier 
for the “Antennae” galaxies. The location of the bend, however, occurs about 1.6 mag fainter in M51. We confront 
the obs erved LF with the model for the evolution of integrated properties of cluster populations of iGieles et alJ 
i2005lJl . which predicts that a truncated cluster initial mass function would result in a bend in, and a double 
power-law behaviour of, the integrated LF. The combination of the large field-of view and the high star cluster 
formation rate of M51 make it possible to detect such a bend in the LF. Hence, we conclude that there exists a 
fundamental upper limit to the mass of star clusters in M51. Assuming a power-law cluster initial mass function 
with exponentional cut-off of the form N dM oc M~^ exp{—M/Me) dM, we hnd that Me = 10® Mq. A direct 
comparison with the LF of the “Antennae” suggests that there Me = 4 x 10® Mq. 


1. Introduction 

There is a relation between the luminosity of the brightest 
star cluster in a galaxy and th e total number of clusters 
llWhitmorell200,3l: lLarsenl^2002^ ■ suggesting that sampling 
statistics is determining the luminosity of the most lumi¬ 
nous cluster. Since the luminosity of clusters is heavily 
dependent on the age, a straightforward translation from 
most luminous to m ost massive is not possible. Recently, 
iHunter et al ] il2nn,'^i showed that the maximum cluster 
mass increases with log(age/yr) in the LMC and SMC, 
which can be i nterpr eted as a size-of-sample effect. Also 
IWeidner et al ] il2nn4 suggest that the maximum cluster 
mass in a galaxy depends on the star formation rate in 
the galaxy, hence the total number of clusters. This sug¬ 
gests that it would be physically possible to form a super 
massive clus ter such as W3 in NG C 7252 with a mass of 
8 X 10^ Mq llMaraston et alJEoOdll in our Milky Way, but 
the chance is just very small. This issue is still heavily 
under debate and is subject of this study. 

The cluster luminosity function (LF) is a powerful tool 
for the study of star cluster populations. In a wide vari- 
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ety of environments the LF can often be well approxi¬ 
mated by a power-law distribution: N dL oc L~°‘dL, 
where the index a is be tween 1.8 and 2.4 (e.g. iLarsenl 
l2nn2t ldr Criis et !m]| 2^. The shape of the LF is related 
to, but not necessarily identical to the cluster initial mass 
function (CIMF). It is important to note that it is hard to 
relate the observed LF directly to the underlying CIMF, 
since the LF contains clusters of different ages. A star 
cluster fades about 5 magnitudes in 1 Gyr in the U-band, 
which makes it hard to estimate the mass without knowing 
the age. 

The LF of clust ers in the “Antennae” galaxies 
I Whitmore et alJll999(). however, is much better approxi¬ 
mated by a double power-law distribution. The bright side 
{Mv ^ —10) has a steeper slope (~ —2.7) than the faint 
side (~ —2). The latter is close to the value found for 
other galaxies. This double power-law nature with a bend, 
was interpreted by the authors as a turn-over in the mass 
function. 

In ICieles et al.l ll2nn5bll we compared a cluster popu¬ 
lation model with various observed luminosity functions 
from the literature. We investigated various possible ways 
of detecting a truncated cluster initial mass function and 
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the possible biases caused by extinction, disruption, vari¬ 
ations in the cluster formation rate, etc. We concluded 
that a truncated CIMF will be observed as a bend in the 
integrated cluster luminosity function. We showed that 
tentat ive hints for a truncation are present in NGC 6946 
ffrom Iba,rsenll2nn^ and M51 Ifrom f^astiau et alJl2nnfijl 
and are clearly not pres ent in the SMC and the LMC 
(from iHunter et al] |2003jl . 

In this work we present a greatly improved LF of clus¬ 
ters in M51, based on recently released deep HST observa¬ 
tions with the Advanced Camera of Surveys (ACS) cover¬ 
ing the entire disk of M51. The great resolution of the ACS 
camera is exploited by selecting clusters based on their 
size. With this we are able to accurately select clusters, 
even when they are as faint as individual bright high mass 
stars. The improved resolution and larger field-of-view 
make it possible to confirm the suggestion of lCieles et alJ 
(1200511 that the LF of M51 is of a double power-law na¬ 
ture. 

We show that the bend in the LF is not necessarily 
related to a corresponding turn-over in the MF, but re¬ 
sults naturally if the CIMF is a power-law distribution 
truncated at the high-mass end. 

In §□ we describe the data, source selection and pho¬ 
tometry. In § 121 we present the LF in the three available 
ACS filters of all extended objects in M51. A comparison 
with the model is done in § 0] and a discussion and the 
conclusions are presented in §[21 


2.3. Photometry 

Aperture photometry was performed on the source list, 
using the IRAF/DAOphot package using a 5 pixel aper¬ 
ture radius and a background annulus with an inner radius 
of 10 pixels and a width of 3 pixels. Aperture correction 
from a 5 pixel aperture to 10 pixels (=0.5”) were mea- 
sure d on artificial sources, gene rated by the BAOl ah pack¬ 
age LarsenlirOOOt iLarsenlBoOdll . A Moffat profile dMoffad 
Il96f ) with power-law index of —1.5 and an effective radius 
(i?eff), which is the radius containing half the cluster light 
in projection, of 3 pc was convolved with the filter depen¬ 
dent point spread function (PSF) of the ACS camera and 
the aperture correction was measured. The PSF we used 
was observationally determined from a crowded star field 
on a drizzled image of the Galactic globular cluster 47 Tuc. 
For each filter a separate PSF was determined. The result¬ 
ing aperture corrections in F^SBW, F555W and F 8 I 4 W 
were —0.16, —0.16 and —0.17 mag, respectively. These val¬ 
ues would be 0.04 lower/higher for sources which are 1 pc 
bigger/smaller. The aperture correc tions between 0 . 5” ari d 
infinity were taken from Table 5 of lSiriau'ni' et all (|2005ll . 
Finally, a correction for Galactic foreground extinction of 
E(B — V) = 0.038 was applied according to Appendix 
B of ISchlegel et al.l l|l99Sll . We did completeness tests by 
adding the same artificial clusters as used for determining 
the aperture correction to the image. A high background 
part of the image was used and the resulting 90% com¬ 
pleteness limits in F435W, F555W and F 8 I 4 W were 23.3, 
23.3 and 23.0 mag respectively. 


2. Source selection and photometry 

2.1. Data 

We used the new HST/ACS (Hubble Heritage) data of 
M51 (NGC 5194) and its companion NGC 5195 in F435W 
(2720 sec.), F555W (1360 sec.), F658N (2720 sec.) and 
F 8 I 4 W (1360 sec.). Six pointings, corresponding to 430 x 
610" (=17.5 X 24.8 kpc), cover the entire disk of M51 plus 
the re gion with NGC 5195. For details on reduction we 
refer to lMutchler et al and the M51 mosaic website 

(http: / / archive.stsci.edu/prepds/m51/). 


2.2. Source selection 


Source selection was don e with the SExtractor package 
dBertin fc Arnout3 Il996t version 2.3.2). A background 
map was made, by computing a mean and standard de¬ 
viation of every section of the image with a user defined 
size. Deviating pixels were iteratively discarded until every 
pixel value of the background was within ±3tT of the mean 
value. Every area of at least three adjacent pixels that ex¬ 
ceeded the background by at least 5cr was called a source. 
For the details of t he background estimation and source se¬ 
lection we refer to lBertin fc Arnoiitj 1 1 996h . The F435W, 
F555W and F 8 I 4 W coordinates of the sources were cross- 
correlated and only sources within a two pixel uncertainty 
were kept. This resulted in a list of 75436 sources. 


2.4. Radius measurements 

To distinguish between stars and clusters, we exploit the 
resolution of the ACS camera (1 pixel = 0.05”) to mea¬ 
sure the radii of all sources detected in F435W, F555W 
and F 8 I 4 W using the ISHAPEroutine within the BAOlab 
package. In summary, analytic profiles with variable effec¬ 
tive radii are convolved with the PSF, and are then fitted 
to each source in the data. The best-fit i?eff was deter¬ 
mined by minimizing the We choose a Moffat profile 
with a power-law index of —1.5. Gomparing the xt of ^ 
point source (= PSF) fit to an extended profile fit, indi¬ 
cated whether a source is resolved. ISHAPE is able to pick 
up sources with a FWHM of 0.2 pixels, which corresponds 
to a ReS= 0. 5 pc at the distyce o f M51 (distance mod¬ 
ulus = 29.62. I^ldmeier et al.lll99^ The source detection 
and radius fits are used to define a source to be a cluster 
when 

1. the source is detected in F435W, F555W and F 8 I 4 W, 

2. the source magnitude is above the 90% completeness 
limit, 

3. the sources is extended, defined as i?efi > 0.5 pc, 

4. the xi of the extended profile fit is lower than that of 
the pure PSF fit. 

These criteria resulted in 5 990 clusters brighter than V = 
23.3 (My = -6.32). 
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Fig. 1. Luminosity function of all extended sources in M51 in F435W, F555W and F 8 I 4 W (left three panels) . The 
solid lines are the double power-law fits (top) and Schechter fits (bottom). The typical bend is indicated as M* for the 
Schechter fits and Mbreak for the double power-laws. Also the slopes are indicated. The dotted line is the extrapolated 
distribution of the faint end. Right: LF of a synthetic cluster population with a CIMF with exponent —2, a mass 
function truncation at Me = 10®Mq and the same detection limit as the data. The top panel shows a double power-law 
fit, where the slope bri ght-ward of the bend is steeper than the underlying mass function. The bottom panel shows 
a fit with a Schechter I Schecht^ll97fi^ function. The thin lines are similar models, but with truncations at 10® Mq 
( higher curve) and lO"* M© (lower curve). 


3. The luminosity function 

Following the source selection as described in § |21 we gen¬ 
erate the LF of all clusters in the three available filters 
(left part of Fig. Pl. The three columns on the left corre¬ 
spond with the F 435 W, F555W and F 8 I 4 W filters. In the 
top panels we fit a double power-law distribution, with a 
variable location of the break and indices and in the bot¬ 
tom panels we fit a Schechter function with variable M* 
and index. The results of the different fits are summarized 
in Tabled Note that the slopes on the faint sides corre¬ 
spond to the typical value of — 2 , found in other galaxies. 
The slopes on the bright side (a 2 ) are much steeper. The 
break appears at higher luminosities in the F 8 I 4 W filter. 
The double power-law break shifts about 0.5 mag to the 
red, going from F555W to F 8 I 4 W. The M* values from 
the Schechter fit are about 0.2 mag brighter in F 8 I 4 W as 
compared to F435W and F555W. 

In iGieles et all il2005lJl we constructed th e LF o f M51 
clusters based on the data of iBastian et al.l ll200,5ll . The 
data was based on two WFPC2 pointings, and the LF 
showed that a double power-law function was better fit 
than a single power-law by a factor of two, when com¬ 
paring xl values. In this case, we obtain xl = 0-78 for 
the double power-law fit in the F555Wha,nd, whereas the 
single power-law fit results in xt = 203. This shows that 
the increase of the field-of-view and depth of the current 
data set contributes considerably to the significance of the 
result. 

The downturn at bright magnitudes and the slopes 
of th e double power-law dist ribution are very similar to 
what IWhitmore et a,lJ (|l 999^1 found for the LF of the 
“Antennae” galaxies. The bend and M*, however, occur 
about 1.6 mag brighter in the “Antennae” LF. We will ex¬ 
plain this behavior of the LF and the difference between 
the “Antennae” and M51 in §0] 


Table 1. Results of the different function fits to the ob¬ 
served LFs of Fig. Hand the simulated LFs of §0]. 



Double power 

-law 

Schechter 

Filter 

«i 

0-2 

Break 

a 

M* 

Observations 

F435W 

1.94 

2.72 

-8.88 

1.81 

-9.87 

F555W 

1.93 

2.75 

-8.93 

1.80 

-9.85 

F814W 

1.93 

3.08 

-9.43 

1.78 

-10.05 


Models with Me = 10® 

Mq 


F439W 

2.01 

2.76 

-8.99 

1.85 

-9.89 

F555W 

1.93 

2.71 

-8.94 

1.79 

-9.94 

F814W 

1.90 

2.70 

-9.60 

1.76 

-10.57 


4. Comparison with a cluster population model 

We have developed an analytical model which can repro¬ 
duce the observable prop erties of a clu s ter pop ulation. The 
model was introduced in iGieles et alJ ll2005all for compar- 
iso n with the age and m ass distribution of M51 and used 
in iGieles et al.l ll2005b|l to compare with LFs in differ¬ 
ent galaxies. The model generates cells which are equally 
spaced in log(age/yr) and log(mass/M 0 ). Weights are as¬ 
signed to each cell, such that the integrated weights as a 
function of mass correspond to a chosen CIMF and forma¬ 
tion rate and in that way correspond to number of clus¬ 
ters. Working with weigths is preferred over generating the 
actual number of clusters, since this can be a very high 
number. Then cells are evolved as if they were clusters, 
where mass loss due to stellar evolution and disruption 
are taken into accou nt by analytical functions derived by 
iLamers et al.l ll2 005ll . 

In iGieles et al.l ll2005bll we showed that a physical limit 
to the maximum cluster mass will cause a bend in the 
LF. For details we refer to that paper. We generated 
several populations with different maximum masses. A 
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constant cluster formation rate between 6 x 10® yr and 
10 Gyr is assumed and a power-law CIMF with an ex¬ 
ponential cut-off N dM oc M~^ exp{—M/Mc)dM, with 
[3 = 2, similar to w hat has been proposed for the CIMF 
of globular clusters ( Burker^^.Smithll200C ). The disrup¬ 
tion parameters from ICieles et all ll2005a 1 are adopted 
(i.e. tdis = 10® (Mi/10^ yr). We note that with 

these disruption parameters all cluster older than 2 Gyr 
are removed from the sample. We applied an extinc¬ 
tion of Av = 0.25 mag to all clust ers, which is close t o 
the average extinction measured bv iBastian et alJ (20^. 
All masses are converted to ACS filter magnitudes, de¬ 


pending on their age, using the GALEV simple stellar 
population models IIAnders fc Fritze-v. Alvenslebenll200.‘lt 
ISchulz et a,lJl2002^ . We fit the same functions to the artifi¬ 
cial LFs as we did to the data. In the right panel of Fig.^ 
we show the result for the F555W filter for Mq = 10® Mq. 
The predicted LF parameters for all filters are summa¬ 
rized in the bottom part of Tabled The bright side shows 
a steeper slope than the faint part, as is the case in the 
data. This is because the clusters with the maximum mass 
fade and therefore have an age-dependent contribution to 
the LF. The brightest cluster will be the youngest and 
going faint-wards, more older clusters will contribute to 
the LF. Faint-ward of My — —9, clusters of all ages con¬ 
tribute. The fact that we find shallower slopes on the faint 
side than the input MF slope is due to the mass depen¬ 
dent disruption. The fact that the bend occurs at a higher 
luminosity in the F814 W filter, in both the model and the 
data, is an important confirmation of our idea that the 
bend is caused by a truncated MF. 

When we make a direct comparison with the 
“Antennae” galaxies, the MF is truncated at a higher mass 
{Me = 4 X 10® AM), since the b end there occurs 1.6 mag 
brighter llWhitmore et al.lll999ll . 


5. Discussion and conclusions 


Does the fact that other galaxies have a single power-law 
LF imply that there is no upper limit to the cluster mass 
there? Probably not. The GIMF has to be sampled well 
enough to reach the critical Mmax) only then a bend will 
show up in the total cluster LF. M51 and the “Antennae” 
are forming enough clusters such that the GIMF is sam¬ 
pled until Mmax, and a bend in the LF is observable. 
The difference between in M51 and the “Antennae” 

galaxies suggest that M^ax is environment dependent. 

These environmental differences might be caused by 
variations in the giant molecu lar cloud (GMG) mass dis¬ 
tribution. < 200, *111 shows that the cloud mass 

function of the “Antennae” galaxies is truncated at higher 
masses than that of M51, whic h in turn is at much highe r 
masses than in the Milky Way dWilliams fc McKeelll 997)1 . 
Other galaxies also seem to have GMG mass distributions 
which are truncated at the high mass end i Posolowskvl 
A truncated GMG mass function might impose a 
physical limit to the maximum star cluster mass, which 


will be observable in galaxies with a high star/cluster for¬ 
mation rates. 
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